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Abstract 

The maximum entropy method has been apphed to investigate the oscillating 
structure in the pp— and pp— elastic scattering differential cross-section at high en- 
ergy and small momentum transfer. Oscillations satisfying quite realistic reliability 
criteria have been found. 
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1 Introduction 



The presence of anomalous structure in the differential cross sections was always indicative 
of some important physical feature typical of the particles and their interactions. Such 
are the dip-bump structure and "fine structure" in the differential cross sections related to 
small and large values of the impact parameter. The first phenomenon is connected with 
absorption at small impact parameters, while the second one to the mesonic (or quark- 
antiquark) "cloud" (or "tail") of the colliding particles. In addition to this, localized 
fluctuations have been observed [|^, |^ in the differential cross-section at fixed y/s. The 
oscillating structure of the diffraction peak in the differential pp— and pp— elastic cross 
sections has been observed for the first time in in the ISR data 0] and later 
in UA4/2 experiment by normalizing the differential cross-section to the smoothly 
varying background in the impact parameter representation. In ^ an attempt was made 
to relate the observed structure near \t\ = 0.1 GeV^ with the variation of the opacity 
in 6— space, probably reflecting the density oscillation in matter. The possible existence 
of oscillations with even smaller periods was discussed by several authors P, ^ §. As 
pointed out in 0, the Coulomb interference region is essential for the clarification of 
the presence of these oscillations in pp— and pp— elastic scattering cross sections. The 
anomalous behaviour observed in the experimental data may be a manifestation of the 
presence of coherent states in the hadronic matter. In ref. ||^ this anomalous behaviour 
was described by a phenomenological formula and interpreted as an unusual phenomenon 
giving rise to new physics in hadron diffraction. 

The smallness of the relevant periods of oscillations indicates the presence of the long- 
range forces. Recently in |^ it was suggested that the long-range forces at t — give 
a significant contribution to the scattering amplitude and may be responsible for the 
structure observed in ^ at small —t ~ 10^'^ GeV^. It has been conjectured in [j^ that the 
oscillations of the differential cross-section observed in at ^/s = 541 GeV are periodic 
in \/i with the period of 2 ■ 10~^ GeV which is compatible with the Auberson-Kinoshita- 
Martin type of oscillations 0. 

In ref. |TDI] we proposed an entirely different method of identifying oscillations in 



the pp— and pp— elastic scattering. Our method was based on the use of overlapping 
bins of local slopes and we indicated some characteristic oscillation periods. It is quite 
obvious that, in order to determine the nature and periods of oscillations, one has first 
of all to increase the reliability of initial information contained in experimental data by 
suppressing the influence of statistical fluctuations. This problem can be settled by means 



of the well known method of maximum entropy |11], used in many branches of physics 



For example, recently it has been applied |T2| to the electron density profile reconstruction 
from multichannel microwave interferometer data at the W7-AS fusion experiment. 



2 The method of maximum entropy 

Consider the mathematical essence of the maximum entropy method. The task is for- 
mulated as follows. Provided there are N measured data points yk with errors Ek and 
no theory available, the question arises of how to choose the most plausible function Yk 
among all possible functions describing this data set. The philosophy is to choose the 
solution which contains the least amount of information in order to avoid false features. 
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Mathematically this means that one has to choose such a function Yk that the functional 
F is minimized: 



F = x' + Ai?. (1) 

Here is the least squares term which provides the consistency of the function Y with 
the data set y, R is a regularization functional which helps to choose the correct function 
containing the least amount of information, and A is the weighting factor between the 
least squares term and the regularization functional. The value of A has to be chosen such 
that the data is neither under- nor overestimated, i.e., should be equal to the number 
of data points N. The is calculated by means of the standard formula 

= E (2) 

k=l 

It is very important to use realistic error values with methods based on regularization 
functional, in order to get realistic results. In the maximum entropy method the regular- 
ization functional is the configurational entropy of the distribution: 

R=-^Pk^og—, (3) 

where 

Pk = > (4) 

ruk - is the fc-th value of the discretised a priori distribution. The basic principle of the 
maximum entropy method is that out of all probability distributions which satisfy given 
constraints, i.e., fits to the data, one should choose the distribution which is closest to the 
given a priori probability distribution m^, and, if this is not specified, one should choose 
the distribution that is closest to the uniform distribution. Here should be pointed out 
that unknowns Y^ are assumed to be independent of each other which means that the 
distribution does not have to be continuous. 



3 Application of the maximum entropy method to 
the differential cross-section 

The question of the role of statistical errors in measured differential cross-section ^ 
for pp— and pp— scattering at fixed s within a broad t-interval is of substantial impor- 
tance, because these errors for each data point may be of the same order of magnitude 
as the "fine" anomalies of ^ observed in experiments 0, ^, 0. The second problem 
is related to ambiguous use of different models and theoretical predictions H, |^ in 
fitting the experimental "curves". From the standpoint of the "observed" oscillations, 
several methods are now being extensively used. One of the methods is based on the- 
oretical description of smooth ^ and subsequent extraction of oscillations (anomalies) 
by subtracting the slowly varying "background" from the experimental differential cross 
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Figure 1: Experimental values of the data at 19.4 GeV ||13[ calculated by means of Eq. 
(jry) (points with error bars), maximum entropy fits with original errors (solid curve; 
X^/ndf = 1.08 and A = 1.5 ■ 10^). 
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For this purpose we minimize the functional (0), where: 



exper 



(5) 



(6) 



/ J \ exper , 

and Ek are experimental errors of 7/^ = ( . We use the parametrization of ^ from 



| 1^ which well describes background to calculate the rrik values needed in (^ 

theor 



da' 



'Fat + Fcf 



(7) 



Fn = CTtot {i + p) exp 



htk + 7 (l - \/l-4Ao 



Here atot, P, b, to and 7 are fitted parameters and Fc is the standard Coulomb amplitude 



which can be calculated from [15]. We investigate the experimental differential pp— and 
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Figure 2: Same as Figure 1 but for the data at 53 GeV 0. Here x^/^df = 1.01 and 
A = 0.23 • 10^ (solid curve). The dashed curve corresponds to the fit with artificial errors 
(three times smaller), x'^/ndf = 0.98 and A = 0.35 • 10^. 



pp— scattering cross sections at different energies covering a broad interval of y/s from 
19.4 GeV |13[ up to 541 GeV [Q. To minimize the functional (|1|), we used the least- 



square method described in [|TB[. The realistic result was obtained when /ndf = 1 was 



achieved. Before using the a priori choice of rrifc in accordance with (0), we set = 1, 
assuming that nothing is known about In this case the use of the maximum entropy 
method just reproduces the experimental behaviour of ^ with slight changes of the slope. 
The results of the calculations with the a priori distribution described above are shown 
in Figures 1-3 for the ratios 
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and the relevant errors. Note that in all cases considered (at 19.4, 53, and 541 GeV) 
y-^j exhibit those non- regularities which are present in l-^j albeit with smaller 
amplitudes. At 53 GeV in spite of the decrease of the amphtude of oscillations, the 
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Figure 3: Same as Figure 1 but for the UA4/2 data 0. Here x^/ndf = 1.04 and A = 
0.35 ■ 10^. 



structure is still visible, especially the local fluctuation at t ~ —0.1 (GeV/c) discussed 
in P]. In contrast to the first two cases, the data on the differential cross-section at 541 
GeV were taken from the region of interference of the Coulomb and nuclear amplitudes 
used to determine the ratio of the real part of the scattering amplitude to the imaginary 
part of the scattering amplitude p{s,t). The maximum entropy method considerably 
reduces the amplitude of oscillations in p{s,t). By means of Fourier analysis we find 
some characteristic periods which generally confirm the calculation of 0, §]. However, 
the value of p [s, t) calculated on the basis of data processed by using the maximum 
entropy method leads to such a drastic decrease of the oscillation amplitude in p (s, t) 
(see Figure 4), that this gives one almost no chance to obtain the value of p{s,t = 0) 
which noticeably differs from the experimental value of 0.135 (compare Figure 3 in |^ 
and Figure 4 of the present paper). 

Finally we made a numerical experiment by artificially reducing three times the experi- 
mental errors of the data of 0. In this case the amplitude of oscillations increased almost 
by an order of magnitude (Figure 1) which stresses the importance of efforts towards 
decreasing experimental errors in real measurements. 
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Figure 4: The ratio p (s, t) as a function of \/t at 541 GeV (bars) calculated directly from 
the experimental data with help 
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The solid line shows the calculations with the differential cross-section corresponding to 
the solid line in Figure 3. 



4 Concluding remarks 

We conclude that in all cases analysed the existing data reveal oscillations satisfying quite 
realistic reliability criteria. Further studies for the small-period structure in t at RHIC 
T7| and LHC (see TOTEM project [^) are highly desirable. To this end good statistics, 



including measurements within the entire diffraction cone with a large number points 
(~ 100) and small errors are necessary. 



The authors are grateful to L. L. Jenkovszky for his constant interest in this work. 
Fruitful discussions with A. V. Snegursky are gratefully acknowledged. 
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